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ABSTRACT 

^ ■ The Australia Telescope Compact Array (ATCA) has been used to survey at 1.4 GHz, 

ly^ ' a small region (< 3 sq deg) overlapping with the 2dF Galaxy Redshift Survey 

l|Colless et al.l200H) . We surveyed with a varying radio sensitivity, ranging from 1 mJy 
Cn • - 20 /xJy (Icr). There are 365 2dFGRS sources with z > 0.001 lying wi thin the sur- 

I veyed region, of which 316 have reliable spectral classification. Following ISa dler et alJ 

(|2002^, we visually classified 176 as AGN or early-type galaxies, and 140 as star- 
forming galaxies. We derived radio flux density measurement or upperlimits for each 
of the 365 2dFGRS sources. The fraction of radio detected 2dFGRS star- forming galax- 
f~i i ies increases from ~ 50% at ^ 0.7 mJy up to ~ 60% at ~ 0.2 mJy. The mean redshift 

Oh' for the fraction of radio detected star-forming galaxies increases with increasing radio 

detection sensitivity, while the mean redshift is fairly constant for the AGN/early-type 
^ • fraction. We found very similar radio detection rates of 2dFGRS galaxies for both the 

^ I AGN/early-type and star-forming components. The radio detection rate increases ap- 

. proximately linearly with respect to the rate of increase in radio detection sensitivity. 

We derived the radio luminosity function for our sample and it was found to be consis- 
tent with that of ISadler et al.l l|2002j) . We have also compared the total flux densities 
of NVSS sources common to our survey, and we discuss strategies for a large-scale 
' radio survey of the 2dFGRS sample. 
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1 INTRODUCTION 

The decimetric synchrotron emission from galaxies is be- 
lieved to arise from relativistic electrons accelerated by 
shock waves either in the interstellar medium, or in struc- 
tures associated with an active galactic nucleus (AGN). In 
the former case, there is evidence that the decimetric lumi- 
nosity is proportional to the rate of massive star formation, 
the link arising because supernovas are p erhaps the dom- 
inant source of accelerating shock waves jBiermannI 119761: 
[Cond on 199^. In the later case, the decimetric luminosity 
of an AGN may be related to the physical properties of the 
underlying black hole, an d the rate and mode of mass accre- 
tion by that black hole llFranceschini. Vercellone fc FabianI 
Il998l: lLaoJl2000l: lHcll200a) ^ 

Observations of the total decimetric luminosity of a 
galaxy thus allow an estimate of the extinction-free star- 
formation rate in galaxies where AGN emission can be con- 



sidered negligible. Conversely, where the star formation rate 
is relatively low they provide information on the character 
of activity excited by a nuclear black hole. A few galax- 
ies emit significant decimetri c from combined nuclear star- 
formation and AGN activity (iHill et al.ll200ltlVeilleuxl200ll : 
iNaear et aL.,,2002'1 . 

The knowledge provided by radio flux density measure- 
ments is greatly enhanced when optical spectroscopic data 
are also available. Optical spectra reveal the galactic red- 
shift, allowing the calculation of luminosities and hence es- 
timates of the cosmic evolution of the rates of star-formation 
and black hole mass accretion. They also reveal the age and 
other properties of the most optically luminous components 
of the galactic stellar population and, in some cases, the 
state of excitement of a substantial fraction of the inter- 
stellar medium of the galaxy. Investigations of large galaxy 
samples having combined optical and radio data offer oppor- 
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tunities to address several outstanding astrophysical prob- 
lems. 

The 2d F Galaxy Redshift Survey (2dFGRS, 

IColless et al.l 1200 has measured the optical spectra 
of over 220,000 galaxies photometically selected to have 
an extinction-corrected magnitude brighter than bj=19.45. 
The survey comprises strips at the equator, and near the 
South Galactic Pole (SGP). The SGP strip densely covers 
~ 80° X 12° (12'^40'" < a < 03''40™, -36.5° <S < -24.5°), 
alo ng with 99 randomly placed 2° (diameter) fields. 

ISadler et al.l (120021) studied the relatively bright radio 
sources among the 2dFRGS objects by cross -matching the 
1.4 G Hz NRAO VLA Sky Survey (NVSS, ICondon et all 
Il998i') with the first ~ 20% of the 2dFGRS. Approximately 
1.5% of the 2dFGRS galaxies (~ 5% NVSS sources) have 
radio counterparts at the NVSS radio flux density detection 
limit of 2.8 mjy. Presumably, the majority of unmatched 
(95%) NVSS 1.4 GHz sources are AGNs whose optical 
counterparts are too faint to be detected at 6j=19.45 (c.f. 
ICondonl ligii'). Some 60% of the identifications are AGNs 
(radio galaxies and some Seyferts) and the remainder 40% 
comprised of star-forming objects. ISadler et al.l i2002l) note 
that complete cross- matching of the 2dFGRS and NVSS will 
provide approximately 4,000 radio source spectra, a sample 
large enough to measure radio galaxy evolution to z = 0.35 
and to lo cate the most lumin ous star-forming galaxies to 
z = 0.2. ICondon et al.l ll2002ll have reported evidence at 
the level of two standard deviations for evolution in the 
rate of star formation density between the mean redshifts, 
(z) ~ 0.02 and 0.06, of t heir NVSS/UGC sa mple and the 
NVSS/2dFGRS sample of ISadler et all (|2002h . respectively. 

There is considerable interest in seeking radio iden- 
tifications of 2dFRGS objects below the NVSS fiux den- 
sity limit. Spectroscopy of the optical counterparts of ra- 
dio sources found i n sub-mjy survey s (e.g. Bonn ct^^]^93 
iGcorgakakis et aL| |l9 99|: Gruppioni. Mignoli fc Zamoranil 
[19998 : .Prandon^tam2001 ) suggests that the faint radio 
population will be a mixture of star-forming galaxies and 
AGNs, as it is for the mJy radio identifications. Interestingly, 
1.4 GHz radio source counts rise sharply towards the Eu- 
clidean rate near 1 mJy (e.g. lCondoJl988l) . a phenomenon 
ascribed to the rising pro portion of star-fo rming galaxies 
in the sub-i njy population (|Wi ndhorstj ll984l : ICondon 1988; 
iBenn et al] ll993': 'Hop kins et alJll998ll . In view of the tight 
correlation between radio and far infra-red (FIR) lumi- 
nosity among the star-forming population, these galaxies 
likely correspond to the population of lower star formation 
rate "IRAS galaxies" seen at apparently brighter magni- 
tudes. Studies of the sub-mJy counterparts of the 2dFGRS 
objects will thus encompass star-forming objects down to 
~ 10^^'^ WHz^^ out to redshifts of z ~ 0.1, where evolu- 
tion appears to have occurred, as well as AGNs of moderate 
radio luminosity. 

The SGP part of 2dFGRS covers over 900 square de- 
grees and would require a very large allocation of radio 
telescope time to undertake a co-extensive sub-mJy survey. 
Acknowledging this, we have undertaken a pilot study of 
a selected region of the 2dFGRS (equivalent area of one 
2dFRGS-field), employing the Australian Telescope Com- 
pact Array (ATCA) at 1.4 GHz to explore probable out- 
comes and optimal observing strategies for a more extensive 
survey. 



Section §2 describes the survey field selection, obser- 
vations and data reduction and §3 the procedures used to 
derive radio fiux measurements and radio upper limits of 
the 2dFGRS objects. §4 presents the radio luminosity func- 
tion derived from our sample and comparison with that of 
ISadler et al.lll2002l) . In §5 we explore the radio detection rate 
of 2dFGRS sources as a function of survey depth and in §6 
we consider the optimisation of a full-scale sub-mJy survey. 

Throughout the paper we have used f^M ~ 1, Ho ~ 
70kms~^Mpc~^ , and adopted for all sources, a radio spectral 
index a = -0.7 {S xu'^). 



2 SURVEY STRATEGY, FIELD SELECTION 

One aim of this pilot survey is to investigate a survey strat- 
egy which yields optimal benefit from the ATCA observ- 
ing time that would be invested in a deep, large-scale radio 
follow-up of the 2dFGRS. The basic trade-off lies between 
the steeply rising surface density of sources at fainter survey 
limits, versus the quadratic dependence of survey sensitiv- 
ity on integration time. This implies a strategy in which a 
small area is surveyed very deeply, and a progressively larger 
area is surveyed to a progressively poorer sensitivity. Given 
the breadth of the total radio luminosity functions for both 
star-forming and AGN galaxies, and the weak dependence 
of the bi-variate radio luminosity function on optical lumi- 
nosity, it would be hard to argue that any flux density range 
is to be preferred over any other, suggesting a strategy in 
which equal numbers of sources are detected to a progres- 
sively changing flux density limit. 

Table 1 presents the results of calculations designed to 
illustrate this strategy. The Table uses the published sen- 
sitivity of the ATCA at 1.4 GHz using two 128 MHz IF 
channels, namely that a point source of Sum mJy is de- 
tected at the held centre with 5a sensitivity in a time of 
{2.24:/ Siim)^ minutes, provided that full {u,v) coverage is 
obtained (ATCA on-line sensitivity calculator^). The Ta- 
ble relies on integral source den sities read from Fig. 2(b) in 
iBecker. White fc Helfaiidl ([199^). According to the Table, a 
single pointing observed for 5,000 m (83 h) would reveal just 
over 3,000 sources per square degree, brighter than 30 fj,Jy, 
and approximately the same number of sources would be de- 
tected to a 5(7 limit of 1 mJy if an area of 34 pointings was 
surveyed for a total of 170 m (3 h). These values are ideal 
estimate, excluding observation overheads. In practice, there 
are several other factors that must be considered in design- 
ing the survey. The limits given in the table refer to the most 
sensitive part of the primary beam, and can be approached 
to within a factor of about y/2 over a large area only if 
mosaicing techniques are exploited. Observations must be 
staged to provide adequate {u, v) coverage, and bright radio 
sources must be avoided to achieve high dynamic range. 

Awarded 8 x 12 hours of observation time for the pilot 
survey, we chose a mosaic of 16 pointings in square tessel- 
lation with a pitch of 17 arcmin. Embedded within the 16 
pointing mosaic are sub-mosaics of 8, 4 and 2 pointings. 
Fig. 1 illustrates our mosaic pattern, where denser filling 
indicates a longer integration time at that pointing. This 
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Table 1. Observing strategy. 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Sensitivity 


Icr 


time 


N > S 


Relative 


Relative 


(MJy) 




(min) 


(deg-2) 


area 


time 


>1000 


200 


5 


90 


34 


1 


>316 


63 


50 


275 


11 


3 


>100 


20 


500 


900 


3.4 


10 


> 31 


6 


5000 


3000 


1 


30 



(1) Point source 5(t sensitivity. 

(2) 1(T sensitivity. 

(3) Time to reach Icr. 

(4) Surface density of 1.4 GHz sources brighter than 5cr limit. 

(5) Relative areas for equal source numbers. 

(6) Relative observing time for specified sensitivity and source 
numbers. 



strategy is not a precise match to the data in Table 1, but 
it sufflces to illustrate and explore it. 

There are several requirements to consider in selecting 
the field. Firstly, it should he away from strong sources (5* > 
100 mjy), and be devoid of sources > 50 mjy in the deepest 
region in order to achieve the desired sensitivity. The NVSS 
has been used to identify regions relatively free of strong 
sources. Secondly, it must lie within 2dFGRS and, for testing 
purposes, have at lea st 50% redshift com pleteness in the 
100k 2dFGRS release ijColless et al.ll200ll ). Finally, we note 
that the 2dFGRS SGP strip passes through the zenith of 
the ATCA (latitude -30°18'52"). Since mosaicing near the 
zenith results in large slew times for the alt-az mount, we 
choose to observe as far away from 5 ~ —30° as possible. 

Using these criteria we place the survey centre at a = 
23''37"35%5 = -27°52'45" (J2000), coincident with the 
2dF fields SGP192-SGP194 & SGP270-SGP272. Fig. 2 dis- 
plays two adjoined NVSS mosaics (I2344M28 & I2328M28) 
showing the neighborhood of the target area. Annotated on 
this image are (i) large circles scaled to represent the 16 
pointings of the 33 arcmin FWHM of the 1.4 GHz ATCA 
primary beam, (ii) small circles representing NVSS sources 
with 50 mJy < S < 100 mJy, (iii) squares representing NVSS 
sources with 100 mJy < 5* < 150 mJy and (iv) numbered 
squares representing NVSS sources with S > 150 mJy. All 
NVSS sources with 5 > 50 mJy lying within a 2° radius of 
the ATCA mosaic centre are marked in the image. 



2.1 Observations and data reduction 

Observations were performed over 8 days at the end of May 
2001 with the ATCA in the 6F configuration. We used 2 x 
128 MHz IFs centred at 1344 MHz and 1432 MHz. Denoting 
the 16, 8, 4 and 2 pointing sub-mosaics as fields A, B, C and 
D, respectively, Table 2 list the average total integration 
time for each of these fields. To simplify observing programs, 
A-fields were observed only during a single 1 x 12 h session, 
B-fields only in 2 x 12 h sessions, C-fields in 4 x 12 h sessions 
and D-fields observed in all 8 sessions. 

To ensure good image quality, we need adequate {u, v) 
coverage. The 6 km array rotates to an independent (u, v) 



Figure 1. Illustration of mosaic pattern and survey depth, where 
each concentric circle groups represents individual pointing cen- 
tres. In accordance with the denotation in §2.1, from lesser to 
denser filling circles indicating increase of total integration times 
corresponding to fields denoted A-fields up to D-fields. 
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Figure 2. NVSS image. Large circles shows the pointings 
(FWHM) surveyed with the ATCA. We have marked all NVSS 
sources within a 2° radius of the ATCA mosaic centre hav- 
ing 50 mjy < 5 < 100 mjy with small circles, and 100 mjy 
< 5 < 150 mjy sources with small squares. Sources with S > 
150 mjy are explicitly indicated. The greyscale spans -10 to 
40 mjy beam~^. 



cell in about 80 s and ideally we would take two samples 
within this time. With 10 s integrations and the relevant slew 
times, this can be achieved only when mosaicing at most 4 
pointings, especially since slew times increase for fields that 
transit near the zenith. Since we have 16 pointings in the A- 
fields, all obtained on a single day, they are under sampled. 
However the other fields are properly sampled since we have 
multiple days of observations to fill the (m, v) plane. 

For primary calibration we used the ATCA standard 
PKS B1934-638, adopting 5'i384 mhz = 14.9 Jy |Baars et all 
1977). For secondary cali bration we h ave used B2331-240, 
adopting Si. 5 ghz = 1 Jy (J^^oJioO^)- The mosaic cadence 
is broken for secondary calibration every 40-50 m, scanning 
the calibrator for 2 m. 
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Table 2. Observation summary, with field-labels defined in §2.1, 
and illustrated in Fig. 1. 



Field 


T* (hr) 


% flagged 


(/xjy/beam) 


A 


0.45 


8 


105 {0.9K) 


B 


1.62 


5 


55 {0.5K) 


C 


6.32 


5 


28 {0.2K) 


D 


25.30 


5 


14 (O.IK) 



* Averaged integration time. 

f Theoretical rms noise level and brightness temperature sensi- 
tivity. 

2.1.1 Calibration and flagging 

The 8 days of observations and each of the two IFs were 
reduced separately following the standard MIRIAD data 
reduction procedures: 

(i) Primary flux calibration: The primary calibrator {u, v) 
data set were flagged to remove bad {u, v) samples and the 
gain, bandpass (mfcal) and leakage (gpcal) calibration ta- 
bles were calculated using a 10 s solution interval. 

(ii) Secondary phase calibration: As the target field was 
near the zenith there is limited coverage of parallactic angle 
for the secondary calibrator. Thus, the bandpass and leakage 
calibration tables from the primary calibrator are used (gp- 
COPy). The secondary calibrator {u,v) data set were then 
flagged and the gain calibration determined (gpcal) using 
a 10 s solution interval. 

(iii) Program source: Calibration tables are copied to each 
of the program sources (pointings) from the secondary cali- 
brator (gpcopy) and each of the (u, v) data sets flagged. 

For our data flagging, we used the task tvclip devel- 
oped bv lPrandoni et al] ll2000l) (see their §4.1). This allows 
the user to recursively flag a {u, v) data set by specifying 
(i) a specific quantity from which a running median of the 
data sample is determined (e.g. amplitude of a polarisation) 
and (ii) a clip level, so that any {u, v) data point exceeding 
the running median by more than the factor specified by 
the clip level will be flagged. We have found that deleting 
(m, v) points having an amplitude greater than 10 times the 
running median removes typical interference. Where lower 
level interference persists we repeat the flagging procedure 
and remove all samples greater than 6 times the running 
median. Table 2 lists the averaged (all days and both IFs) 
percentage of data flagged for the field types defined in §2.1. 
The overall total amount of data flagged is 4.5%. 



2.1.2 Imaging and self-calibration 

Once the data are calibrated and interference is ex- 
cised, we proceed to masaicing. We have adopted a joint 
deconvolution appr oach with a Steer-Dewdney -Ito (SDI) 
CLEAN algorithm (|Steer. Dewdnev fc Itolll984h using the 
task MOSSDI 2 . This task is a version of mossdi modified by 
iReean et all J200ll) so that the CLEAN-cutoff level is ex- 
pressed in terms of the locally determined image RMS in- 
stead of a fixed value in Jy beam~^. This feature was very 
important in view of the large variation in the total inte- 



gration time per pointing across the mosaiced region. We 
followed the basic MIRIAD joint-deconvolution procedures: 

(i) Apply INVERT to the {u, v) data of all pointings, both 
IFs (invoking options mfs to deal with bandwidth smearing) 
and all 8 days, with uniform weighting, 2 arcsec cell size, 
to generate a dirty mosaic and dirty beam cube (one dirty 
beam image per pointing). 

(ii) Use MOSSDI2 to generate a joint-deconvolution 
CLEAN-component model of the dirty mosaic, CLEANing 
down to a specified CLEAN-cutoff level. 

(iii) Use restor to subtract the CLEAN component 
model, convolved with the appropriate dirty beam, from the 
dirty mosaic and then restore with an elliptical Gaussian 
profile fitted from the dirty beam. 

(iv) Lastly, use mossen to generate a gain image which 
then allows primary beam correction over the entire mosaic. 

There are two artefacts evident in the mosaic. Firstly, 
weak positive and negative rings appear around strong 
sources, presumably reflecting amplitude calibration errors. 
The effect is largest for bright sources lying on the flank 
of the primary beam and is probably caused by pointing 
errors. There are also radial artifacts appearing as symmet- 
ric and asymmetric radial structures emerging from strong 
sources. These can be removed by phase self calibration. 
Comparison of our mosaics before and after the self cali- 
bration procedures, we found no spurious sources above 3a 
level, and source flux densities varies less than ~ 5%, much 
smaller than the expected fitting uncertainty. The following 
summarises our imaging steps: 

(i) Without restricting the CLEAN region deconvolve 
with MOSSDI2, CLEANing down to a CLEAN-cutoff level 
of 10(7 and generate a preliminary mosaic. 

(ii) Use automatic source fitting of the preliminary mosaic 
to generate a list of CLEAN boxes over all fitted sources with 
peak flux density greater than 10 times the local image noise 
rms. 

(iii) Generate a CLEAN component model for self- 
calibration using MOSSDI2, CLEANing over the boxes de- 
flned in (ii). Invoke the option positive to constrain the 
CLEAN components to be positive valued. 

(iv) With the CLEAN components deflned in (ii) in- 
vestigate image quality improvement using various self- 
calibration solution intervals. Radial artifacts were com- 
pletely removed for solution intervals lying between 10-30 s. 
We adopted 30 s to enhance the signal-to-noise ratio per 
self-calibration solution interval. 

(v) With SELFCAL and a 30 s solution interval apply phase 
self-calibration to each pointing, each IF, and each day sep- 
arately. 

(vi) Using self-calibrated (u, v) data, deconvolve with 
MOSSDI2 CLEANing down to 5a over CLEAN boxes deflned 
in (ii). Use the new mosaic to define a new set of CLEAN 
boxes for all fitted sources with peak flux density greater 
than 5 times the local image noise rms. Finally generate a 
deep-CLEANed mosaic with MOSSDI2, CLEANing down to 
2a over the new CLEAN boxes. 

Displayed in Fig. 3 is the deep-CLEANed mosaic with 
an enlarged central portion. The contours represent la local 
image noise levels calculated in the manner described below. 
The 0.2 mjy noise contour in Fig. 3 corresponds closely to 
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Figure 3. ATCA mosaic. Tiie contours are drawn at tlie Icr noise 
level for 0.025, 0.05, 0.1, 0.2, 0.4 and 1 mjy. The greyscale is set to 
-1.5 to 4.5 mJy beam~^ for full the mosaic and -0.15 to 0.45 mJy 
beam~^ for the inset. 



the outline of the 16-pointing mosaic region illustrated in 
Fig. 2. The synthesised beam size of the deep-CLEANed 
mosaic is 15.5" x 8.0" at position angle 1.20°. 

2.2 Local image noise estimate 

It is important to estimate the local RMS noise for the im- 
age used in the source-fitting phase of our study. We could 
estimate the noise using restor to form a residual image of 
the deep-CLEANed mosaic, but the estimate would exhibit 
irregular spatial variations. Instead we have generated a the- 
oretical sensitivity image using mossen, and checked the 
predictions using noise estimates based on 1000 randomly 
selected patches each ~ 3x the synthesised beam area. The 
tight correlation between the prediction and the measured 
values shown in Fig. 4 is fit by a = lO"""*^ {n)°"'' with ~ 
20% uncertainty. This equation is used to scale the theoreti- 
cal sensitivity image to estimate the local image RMS noise. 
Table 2 lists the theoretical RMS noise of the four field types 
defined in §2.1: the raw mosaic image RMS noise agrees with 
the predictions to within ± 20%. 



3 RADIO DETECTION OF 2DFGRS 
GALAXIES 

In measuring the radio properties of 2dFGRS galaxies in our 
pilot survey, we restricted our radio image to a lowest RMS 
sensitivity of 1 mJy, the outermost contour in Fig. 3. This 
is necessary as to limit the uncertainty in the primary beam 
corrections of our radio image. The total area to a 1 mJy 
detection limit is ~ 3 deg^. Fig. 5 shows the variation of 
survey area as a function of la sensitivity. 

Based on the final 2dFGRS data release, a total of 392 
2dFGRS objects he within the 1 mjy-region of the ATCA 
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Figure 4. Plot of residual image rms against theoretical sensitiv- 
ity image mean, with a best-fitting line (solid) and approximated 
range of uncertainty (dashed-lines). 
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Figure 5. Survey area as a function of Icr RMS noise level. The 
solid line is a cubic spline interpolation of the data. 



mosaic. Of these, 27 have reliable (Q > 3) low-redshift spec- 
tra (2: < 0.001), they are stars and have been excluded 
from further analysis. The spectra of the remaining 365 ob- 
jects have been visually classified by author BC and cross- 
checked by author E MS. We adopted the visual classifica- 
tion scheme used by |Sa^«^t al. (200 3) a subset of which 
were also classified bv lJackson fc LondishI (|200(J) using di- 
agnostic emission-line ratios. The spectra were classified as 
AGN galaxies with (1) ahs: absorption line spectra, (2) em: 
emission-line spectra with dominant forbidden lines (e.g. 
[QUI]), or (3) a+e: combined absorption and emission line 
spectra, or as star-forming galaxies with (4) SF: emission- 
line spectra with dominant recombination lines (e.g. Hq), 
or as un-classifiable objects (5) ???: objects with spectra 
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Figure 7. The differential matciies of radio-2dFGRS sources as 
a function of radio-optical offset, shown with Poisson errorbars 
(^/n). Also shown are averaged Monte-Carlo test results and the 
solid line shown the expected chance coincidence level, assuming 
radio sources are uniformly distributed. 



Figure 6. Examples of 2dFGRS spectra, classified from the top 
as (1) abs, (2) em, (3) a+e, and (4) SF. See §3 for details. 



which are too noisy or otherwise un-classifiable. Note that 
due to aperture effects of fiber spectroscopy, a spectra classi- 
fied as a bs or a+e may dep end on the redshift of the galaxy 
(e.g. see lSadler et alj|200l their Fig.4 and §3.4). Fig. 6 ex- 
hibits sample 2dFGRS spectra in classes (l)-(4) from top to 
bottom. 

We aim to determine the radio flux density or an upper 
limit at each 2dFGRS source position. For the brighter radio 
sources, this can be done by cross-matching a radio source 
catalogue to the 2dFGRS catalogue. However, since the op- 
tical source positions are known we can explore potential 
low signal-to-noise (S/N) detections by examining the sta- 
tistical outcomes of source fitting at the known optical posi- 
tions. Our strategy for cross-identification when the optical 
positions are known employs the AIPS procedure VSAD to 
catalogue radio sources with peak fiux density > 5 times the 
local image noise rms, and the MIRIAD procedure imfit to 
examine the statistical properties of pixels near the optical 
sources. Estimated uncertainty in the flux density measure- 
ments by the two source fitting programs have been derived 
using Monte-Carlo simulations, the results are presented in 
Appendix I. 

We deploy vSAD for multiple Gaussian fitting within 
regions of ~ 3 beam areas, centred in the previously deter- 
mined CLEAN boxes and adopting a peak fiux density cutoff 
of 3 times the local RMS noise determined as explained in 
§2.2. The resulting 3a peak flux density list is filtered to 
produce a radio source list with peak flux densities exceed- 
ing the local 5a level and inspected visually to reject obvious 
spurious fits. This radio source catalogue, containing 326 sin- 
gle components and 2 paired sources, is then cross-matched 
with the 365 2dFGRS positions. To determine the optimal 
search radius for cross-matches, we consider Fig. 7 which 
shows the number of matches ("o" with \/7^ Poisson er- 



rorbars) inside consecutive concentric annuli centred at the 
optical position. Monte-Carlo techniques were used to deter- 
mine the level of chance coincidence in the case of variable 
survey sensitivity. The Monte-Carlo results ("*" -symbols) 
agreed well with the expected chance coincidence (dotted- 
line in Fig. 7) calculated for a uniform radio source density. 
Fig. 7 shows that beyond a radio-optical offset of about 4- 
5 arcsec, the rate of radio-2dFGRS cross-matches overlaps 
with that expected from chance coincidence. Matches with 
offsets up to 5 arcsec are accepted as automatic identifica- 
tions, but we visually inspect matches up to 15 arcsec offset 
to uncover special cases. 

Of the 365 2dFGRS positions, we found 23 identifica- 
tions with radio-optical offsets smaller than 5 arcsec, and 3 
positions with 5-10 arcsec offsets which were all accepted 
after visual inspection. The 2 matches with offsets 10-15 
arcsec were rejected as identifications. 

For each 2dFGRS position without a radio match, the 
MIRIAD task IMFIT was used to attempt a least-squares fit 
of the radio point spread function at the optical position. 
Although the fitting algorithm is understandably unstable, 
fits converged for 65 cases out of the remaining 339 optical 
positions, with a fitted peak flux density greater than the 
estimated local RMS noise (Icr). Of these, 9 having a peak 
flux density > 3a are reported as flux density measurements. 
The 56 with fits below 3cr are characterised as upper limits. 
For the remaining 274 optical positions with unsuccessful 
IMFIT results (divergent fits or fltted flux density < la) we 
report a flux density of zero, and specify the estimated local 
noise as the radio upper limit. 

We developed a glish-scnpt to generate image overlays 
using methods in the viEWER-tool of AIPS++. Red 2"^^ 
epoch Digitized Sky Survey (DSSII-red) images'^ were over- 
laid with radio contours at each of the 365 2dFGRS objects. 

^ DSS images obtained from the Canadian astronomy data cen- 
tre, at: http://cadcwww.dao.nrc.ca/dss/ 
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DSSIl-red image with ATCA radio cantourB 




H3^38'^5E" 50° 4a" 

Right Ascension (JSOOO) 
^ I I^SSIT-reri imj^ci; ^ with ^TCfl i-nriio r^ontoiir^ 




: I I I I I ^ I I I 1 1 1 I 1 I I I 

£ai'34"6° 34™4'' 34"3'' 34'^0= aS^SB" SS^Se" 

Right Ascension (J2000) 

Figure 8. DSSII-red image overlaid witli (top) radio contours 
at 2, 3, 5, and 7 Xa = ATA fJ,3j) and (bottom) radio contours at 
2, 6, 11, 17, and 24 Xa = 358 /iJy). The dashed outUne marks 
the region used to determine the radio flux density and the cross 
marks the 2dFGRS position. 



These images led us to accept all 3 matches with radio- 
optical offsets between 5-10 arcsec, and also revealed two 
objects whose radio emission was extended and required 
pixel summation to determine the integrated flux densities. 
DSSII-red images of these objects are displayed in Fig. 8, 
overlaid with radio contours. The dashed-outline in Fig. 8 
shows the region used to determine the radio flux density 
and the cross marks the corresponding 2dFGRS object po- 
sition. 



3.1 ATCA-2dFGRS data-table 

We publish here the radio flux density measurement or 
upper-hmits of the 365 2dFGRS objects within our ATCA 



radio survey region. Table 3 shows page 1 of the full data 

table. This catalogue is sorted on right ascension, and the 

structure of the table is as follows: 

Columns 1 & 2: 2dFGRS position (J2000); 

Column 3: 2dFGRS object name; 

Column 4: 2dFGRS redshift; 

Column 5: final Bj magnitude used in 2dFGRS; 

Column 6: 2dFGRS redshift measurement quality; 

Column 7: visual spectral classiflcaton; 

Column 8: fltting method used to determine the radio flux; 

(1) vSAD-5cr match; 

(2) IMFIT > 3ct match; 

(3) IMFIT 1 — 3 (J match; 

(4) local noise upper limit; 

(5) pixel sum; 

Columns 9 & 10: vsad-5(t fitted radio position; 

Columns 11: estimated VSAD fitted positional uncertainty 

in arcsec, see Appendix I; 

Columns 12: radio-optical angular separation; 

Columns 13 & 14: peak & integrated radio flux density 

determined using the corresponding fitting method; 

Columns 15 & 16: percentage & actual radio flux density 

uncertainty^, see Appendix I; 

Column 11: estimated local noise, see §2.2; 

Column 18: radio signal-to-noise ratio of the detection. 



3.2 Summary and comparison of results 

Of the 365 2dFGRS objects with z > 0.001, 316 have reliable 
spectral classification, with 176 and 140 classified as AGNs 
and star-forming galaxies, respectively. The remaining 49 
objects were unclassified, having either a low quality fiag (Q 
< 3) or noisy spectra. 

Of the 35/365 (9.6%) radio detections (measured ra- 
dio fiux density exceeding 3cr) of these galaxies, 12 were 
classified as AGNs, 20 as star-forming galaxies, and 3 were 
unclassified. Table 4 presents a summary of the radio detec- 
tion rate of AGNs and star-forming galaxies in 2dFGRS for 
th e current survey an d those identified with NVSS sources 
bv lSadler et alj l|20o3). It is diflicult to compare the NVSS 
and ATCA values in view of the varying sensitivity in the 
ATCA survey, but they do reveal a rise in the detection rate 
with deeper radio sensitivity, especially for the star-forming 
population. 

Fig. 9 shows the redshift and radio luminosity distri- 
butions of the radio-identified 2dFGRS sources associated 
with AGN and star-forming populations. To facilitate com- 
parison, each distribution were normalised by the total of 
the distribution. Note that the narrow range of radio lumi- 
nosities found in the ATCA sample is a selection effect owing 
to the combination of optical-radio survey cutoffs and the 
relatively small surveyed area. The sub-mjy ATCA survey 
detects an increasing proportion of lower-luminosity AGNs 
lying over a redshift range similar to that of sources de- 
tected above 2.7 mjy, as well as an increasing proportion 
of higher-redshift star-forming galaxies having a mean lu- 
minosity slightly lower than the NVSS sample. This implies 

When the local noise exceeds the fitting uncertainty, the larger 
value was used. 
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Table 4. ATCA-2dFGRS summary. 





2dFGRS 


ATCA 
detection 


NVSS 
detection 


z-unreliable 


49 


3 (6.1%) 




AGN 


176 


12 (6.8%) 


~ 1.8% 


SF 


140 


20 (14%) 


~ 1.1% 


Total 


365 


35 (9.6%) 


~ 1.5% 


AGN:Zn,^d (sd) 




0.14 (0.046) 


0.14 (0.16) 


SF:Zmed (sd) 




0.086 (0.037) 


0.046 (0.040) 


AGN;logio(Li.4 


Gfl"z)med (sd) 


22.22 (0.69) 


23.58 (0.81) 


SF:logio(Li.4 GHz)med (sd) 


21.93 (0.47) 


22.38 (0.58) 



Redshift distribution Radio iLininosity distribution 




■ — I ' ' ' I ' ' ' ' I I ' ' ^ ' L' 1 ' ^ 
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'°0,„|L,,„H.[™* 11 



Figure 9. Redshift (left-column) and radio luminosity (right- 
column) distributions for AGNs (top-row) and star-forming galax- 
ies (bottom-row) of 2dFGRS sources identified with this current 
radio survey (dashe d-line-histogram) a nd NVSS sources (solid- 
line-histogram) from lSadler et all l2002^ . The two vertical dotted- 
lines in each plot shows the median of the distributions, the actual 
values are listed in Table 4. 



that a sub-mjy radio survey of the 2dFGRS will better re- 
veal the population in the faint end of the AGN and star- 
forming galaxy radio luminosity function, and enable more 
reliable determination of the evolution of the star-forming 
population to z ~ 0.1 (radio hmited, 50 piJy) and the AGN 
population to z ~ 0.3 (optical limited). 



source, the maximum detectable redshift range [zmin, Zmax]- 
Adopting the Einstein-de-Sitter model, the maximum de- 
tectable co-moving volume (Vmax) is given by: 

("Zmax [1 I 1-)-l/2l2 

Vma.^l&Ajj^f j L0{z Y , (1) 

where uj{z) is the fraction of the surveyed solid angle at 
redsh ift z. The adopted o ptical survey limit is 14.0 < b,j < 
19.4. (|Golless et alJ (|200lf l'). constant throughout the survey. 
The radio detection limit varies across the radio mosaiced 
region (see Fig. 3), and is set at 3a. A cubic spline interpo- 
lation of the relation shown in Fig. 5 provides an estimate 
of uj{z) for any particular source, and Vmax is obtained via 
numerical integration of cquation^following the method de- 
scribed in §2(d) of Avni & Bahcall (1980). Similarly, we can 
calculate the detectable co-moving volume of each source 
iy ). The ratio V/Vmax of a randomly distributed sample is 
expected to be uniform in the interval [0,1], and the mean 
{(y/Vmax)) should be ~ 0.5. 

Our sample comprises 35 ATCA-2dFGRS sources which 
satisfy the survey limits stated above. Three without reli- 
able redshifts are excluded from the RLF calculations, but 
included through the incompleteness correction. To estimate 
the level of incompleteness of our ATCA-2dFGRS RLF, we 
examined galaxies in the 2dFGRS input catalogue that fall 
within the radio surveyed region but were not observed with 
2dF. Only one of these "missing" 2dFGRS sources was iden- 
tified with a > 3a radio source, implying that this optically 
selected radio-sample is approximately 90% (32/36) com- 
plete. 

The total RLF and contributions by the AGNs and SF 
galaxies are listed in Table 5, also listed are the {V/Vmax) 
values. While the SF galaxies in our sample seems to be uni- 
formly distributed, the V/Vmax values of the AGNs show a 
small but not significant increase and they are non-uniformly 
distributed {xneA\a.n{V /Vmax) ~ 0.617). This might indicate 
some level of clustering or evolution in the AGN fraction 
within our sampled volume. 

Presented in Fig. 10 is the total RLF of our ATCA- 
2dFGRS sample (*-solid) al ong with val ues derived from 
the NVSS-2dFGRS sample of lSadler et all £2002). The RLF 
separated into the populations of AGNs and SF galaxies of 
the two samples is shown in Fig. 11. The values and the plots 
have not been corrected for the 10% incompleteness since the 
ef fect ( log 1 n ( 1 ■ 1 ) ) is small . Our RLF agrees well with that 
of lSadler et alJ |20o3), within measurement uncertainty and 
the luminosity range sampled, due to the smiliar redshift 
range sampled. Although our median redshift for the star- 
forming galaxies are higher, with such a small sample we 
are not yet able to test for any significant differences with 
respect to other samples. 



4 RADIO LUMINOSITY FUNCTION 

In this section we describe the construction of the total radio 
luminosity function (RLF) and the corresponding contribu- 
tions by the AGNs and star-forming (SF) populations for 
our ATCA-2dFG RS sample. We have used the traditional 
1/Vmax method oflSchrnidti lll968fl to derive the RLF, as de- 
scribed in §2 of lCondonI (|l2S3)- According to the combined 
optical and radio survey constraints we calculate, for each 



5 RADIO-2DFGRS DETECTION RATE 

To explore the radio-2dFGRS detection rate at various radio 
detection limits, the radio survey was divided into four zones 
of sensitivity, starting from 25^Jy and doubling in four steps 
upwards. The four zones thus corresponded to 3a radio de- 
tection limits of 75, 150, 300 and 600 /xJy. The number of 
2dFGRS sources and radio detections above the 3a limit are 
counted in each zone, separated into three optical spectral 
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Table 5. ATCA-2dFGRS radio luminosity function. 





All 


galaxies 


SF 


galaxies 




AGNs 


Wio(i*) 


N 


loglo(Pin) 


N 


logio(Pm) 


N 


logio(/oS) 


20.5 
21.3 
22.1 
22.9 


1 
5 
23 
2 


r, -,, + 0.30 

-3.5o«:;? 

-4 62+°-2'* 


1 
5 
14 


9 cn + O.SO 

"^■"^-1.00 

-3 40+°-^2 
■^■^^-0.48 

-3 67+" " 

■^■"'-0.15 


9 
2 


3 nQ+0.13 

-4 62+°-^* 
^•°^-0.60 


23.7 





<-5.49 









<-5.49 


24.5 


1 


5 76+0.30 






1 


5 7fi+0.30 


Total 


32 




20 




12 






0.526 ± 0.051 


0.492 ± 0.065 


0.584 ± 0.083 


* ii.4GHz in 


[WHz 














"=S,o('-,.4GHz[™"^ 1) 

Figure 10. Tota l radio luminosity f unctions for the 2dFGRS- 
NVSS sample of ISadler et all <2002l) (o dashed-line) and our 
2dFGRS-ATCA sample (* solid-line). 



classes: (i) star-forming (SF) galaxies, (ii) AGNs and (iii) 
no redshift (unknown). These counts are listed in Table 6, 
excluding the dete ctions above the 2.8 mjy limit used in 
ISadler et all ll2002l) . 

Our survey region has been specially selected to avoid 
bright radio sources, therefore the radio-2dFGRS detection 
rates calculated from our sample could be biased against 
identifications with bright sources. To overcome this possi- 
ble bias, we have adopted the expected detection rate for 
sources above 2.8 mJy using the NVSS-2dFGRS results of 
(Sadler et al. 2002). Fig. 12 illustrates the proportion of 
2dFGRS sources identified with radio sources as a function 
of radio survey sensitivity. This is calculated from the sum 
of the detection rates at the various radio detection limits 
listed in Table 6 and the corresponding detection rate from 
the NVSS-2dFGRS sample. 

The plot shows an approximate 10-fold increase in the 
detection rate for a 10-fold improvement in radio detec- 
tion sensitivity. This increase in detection rate would re- 
quire a 100-fold increase in observing time. These results 




Figure 11. Radio luminosity functions separated into the star- 
forming (lighter-colour lines) and the AGN (darker-colour lines) 
populations for the NVSS (o dashed-line) and ATCA (* solid-line) 
samples. 



are broadly consistent with the estimates made in §2, us- 
ing our observing strategy, we can obtain a sample evenly 
distributed in numbers of identifications, independent of the 
additional constraint imposed by 2dFGRS optical selection, 
across all radio flux densities. 



6 LARGE SCALE RADIO SURVEY 

Owing to the enormous time investment in acquiring large 
optical-spectroscopic redshift surveys such as that by 2dF 
makes it even more valuable to fully exploit its potentials. 
A large scale radio follow-up survey of the 2dF sources will 
fill a significant gap in the available data on faint radio 
sources. Existing radio-optical spectroscopic studies have 
been based on either (i) bright, large-scale radio surveys 
(FIRST*, NVSS and SUMSS'^) that are cross-correlated 
with catalogued optical data (implying that the identified 
optical galaxies are also bright) or (ii) deep radio sur- 
veys made over small areas (such as: PDS'', HDF^ and 
HDFS*) with pointed optical spectroscopy of selected op- 
tical IDs. Studies of type (i) have given radio-optical data 
on ~ 4000 and ~ 100 nearby galaxies from the UGC-NV SS 
JCondon et alj|2002l) and 2dF-NySS iSadler et al.ll2002D & 



2dF-FIRST jMagliocchetti et al.ll2000D . with z < 0.1 and 
z < 0.3, respectively. Studies of type (ii) (which is very ex- 



Very Larg e Array-Faint Images of the Ra dio Sky at Twenty- 
centimeters iBecker. White fc Helfandlll995D. 

^ Syd ney University Molonglo Sky Survey jBock. Large &: Sadlerl 

^ Phoenix Deep Survey fHopkins et al."l99i 
Hubble Deep Field ( Richards et al. 199"" 



Hubble Deep Field South iNorris et alTBoOll) 
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Table 6. Radio- 2dFGRS identification counts at the corresponding Scr radio limit up to the 2.8 mjy limit used in lSadler et al.l 



Radio 3(T limit Area Number of 2dFGRS galaxies Number of radio identifications 

(mJy) (deg2) Total SF AGNs unknown Total SF AGNs unknown 



0.075 - 


2.8 


0.19 


27 


8 


17 


2 


7 


2 


5 




0.150 - 


2.8 


0.72 


100 


43 


45 


12 


18 


10 


6 


2 


0.300 - 


2.8 


1.30 


163 


73 


74 


16 


14 


10 


4 




0.600 - 


2.8 


1.93 


227 


97 


101 


29 


6 


4 


2 






Radio survey detection limit [mJy] 

Figure 12. Radio-2dFGRS detection rate at various radio detec- 
tion limits. The lines interpolates the detection rate with Pois- 
son errors for: top-panel - the total-detection rate (dotted-line); 
bottom-panel - detection rates of star-forming galaxies (solid- 
line) and the AGNs (dashed-line). Note: (1) the d etection rate 
at 2.8 mJy is from the NVSS-2dFGRS sample of ISadler et alJ 
l2002h . (2) the detection rates above 2.8 mJy are the sum of 
the corresponding NVSS-2dFGRS detection rate and the values 
shown in Table 6, see text for details. 



pensive of telescope time) has to date provided radio-optical 
data on fewer than 300 galaxies, at redshifts up to 2 ~1. 

Fig. 13 displays the radio-detected 2dFGRS objects in 
the current survey and those in 2dF-NVSS. Also shown are 
dashed lines tracing out the UGC-NVSS sample, and a solid 
line delineating the main concentration of UGC-NVSS data 
points (predominately star-forming galaxies) . It is clear that 
mJy radio surveys are only capable of properly sampling the 
full range of radio luminosities for the SF and AGN popu- 
lations up to 2 ~ 0.01 (I/1.4GHZ 1 10^^'^ WHz"^). Whereas, 
even though the ultra-deep radio surveys can provide radio 



luminosity sampling to very high redshifts {z > 1), the prob- 
lems lies in that pencil beam surveys (i.e. small survey area 
of typically a few square arcmin up to a few square degrees) 
are not ideal to sample galaxy evolution. There are inherent 
problems in small area surveys associated with clustering, 
cosmic variance, and the complex selection effects will be 
both difficult to model and correct for. A large area sub-mjy 
radio surveys of 2dFGRS sources with a varying sensitivity 
will link together the existing surveys and provide the re- 
quired sample of the 'normal' population of sources in the 
intermediate redshift ranges up to z ~ 0.03-0.04. An investi- 
gation of possible systematic source flux differences between 
our ATCA survey and those in the NVSS, we present, in 
Appendix II, flux comparison of sources found common in 
the two surveys within our mosaiced region. 

ISadler et al.l i2002l) estimate that approximately 4,000 
radio-optical identifications will be found when the full 2dF- 
GRS is correlated with the NVSS catalogue (Jackson et al. 
2003, in preparation). Using the detection rates measured 
with our pilot study, we can now estimate the survey sen- 
sitivity and the total survey area (i.e. observing time) re- 
quired to provide a set of radio IDs that is complete to a 
given radio luminosity and redshift. Computed in (Table 7) 
is the required observing time and the number of expected 
optical-radio identifications for a survey with three levels of 
sensitivity and survey area. The three adopted 3a sensitiv- 
ity levels form a bridge between existing surveys, the poor- 
est sensitivity level at 0.6 mJy being approximately half of 
the FIRST limit (the FIRST limit is approximately half of 
NVSS limit) while the best sensitivity level of 0.15 mJy is 
about half as good as that of PDS. 

Employing a nested mosaicing technique the survey 
plan outlined in Table 7 will require a total of about 60-70 
X 12 hour observing sessions. The survey will image about 
100 deg^ of sky, making over 1000 radio/optical identifica- 
tions. This number is required to address the important 
problem of constraining the cosmic evolution of the radio 
spectral power density of star-forming galaxies. The UGC- 
NVSS sample has provided a very precise measurement of 
the local {z ~ 0.01) star formation density, presenting a lo- 
cal reference and has hinted at intriguing cosmic evolution of 
this guantitv lCondon et al.l (I2OO2I) . With favored evolution- 
ary model in which the luminosity density rises as (l-l-z)"^^^, 
the predicted increase of star formation rate density is ~ (30 
± 10) % at 2: ~ 0.1. With a sample size of 1000, divided into 
ten luminosity bins for deriving the redshift-dependence of 
the star formation density, producing a measurement uncer- 
tainty of ~ 10%. This would be sufficient in distinguishing 
between various models for the cosmic evolution of the star 
formation density, providing a measure potentially superior 
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Radio flux limits 



0.001 0.01 0.1 1 

Figure 13. Plot of radio luminosity against redshift, for our 
ATCA-2dFGRS (3o-) detections along with the NVSS-2dFGRS 
(~ 5it) sample. The dashed-line traces out the UGC-NVSS (~ 5cr) 
sample, the solid-line shows the main concentration of data points 
(predominately star-forming galaxies). The dotted- lines shows 
various flux density limits. 



to the traditional estimates via optical emission-line mea- 
surements or ultra-violet photometry which are unavoidably 
biased by dust extinction. 
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APPENDIX I: SOURCE FITTING AND 
UNCERTAINTIES 

In measuring the radio flux density associated with 2dF- 
GRS sources, we have used two automatic radio source fit- 
ting packages: (1) the MPS task VSAD, and (2) the miriad 
task IMFIT. We first matched the optical source list to a 
radio catalog generated with vsad, of bright (> 5(t) radio 
sources. For positions without a bright radio counterpart, we 
used IMFIT to determine radio upper-limits by fitting a beam 
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Table 7. Sample observing plan. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


Area 


# pointing 


So- 


ID rate 


time/field 


# of 12hr 


#ID/dcg2 


#ID below 


Total 


(deg2) 




(mJy) 


{%) 


(min) 


observation 




NVSS/deg 


#ID 


98.0 


1221 


0.60 


5 


7.5 


13 


5.12 


3.6 


351 


49.0 


610 


0.30 


10 


30.0 


25 


10.23 


8.7 


426 


24.5 


305 


0.15 


20 


120.0 


51 


20.47 


18.9 


464 



(1) Area surveyed (in square degrees) to specified sensitivity. 

(2) Number of pointings needed to cover required area. 

(3) Radio 3cr (3 X image noise rms) sensitivity limit. 

(4) Percentage of expected optical-radio identification rate, obtained from pilot study. 

(5) Amount of radio telescope time required to reach sensitivity limit per pointing. 

(6) Number of 12 hour observing sessions required to cover area. 

(7) Number of optical-radio identification expected per square degree. 

(8) Number of optical-radio identification expected per square degree inaccessible to NVSS. 

(9) Total number of optical-radio identification dowrn to the sensitivity limit inaccessible to NVSS. 



size elliptical Gaussians at the known optical positions. Fol- 
lowing these techniques, we investigated the source fitting 
uncertainty with Monte-Carlo simulations. 

Beam-size Gaussians were injected at randomly chosen 
positions into the residual image of the deep-CLEANed mo- 
saic. To obtain an estimate of the source fitting accuracy 
as a function of fitted source strength^ (""a), a range of in- 
jected source strengths were used, corresponding to a signal- 
to- noise ratio up to 150 a. We first used VSAD to generate 
a 5a radio source list, following the same procedures as we 
used in our study (see §3). The input position list is then 
matched against the 5a radio source list to within 10 arcsec, 
while fitting with imfit at the known source position for 
those without a match. 

Shown in Fig. 14 is the Monte-Carlo results for the vSAD 
fitted angular positional deviation from the input at various 
fitted source strengths. The solid-curve is a quadratic least- 
square fit to the logarithmic one-standard deviation binned 
values of the results. This curve bounds about 90% of the 
data points. For instance, for about 90% of the 5a VSAD 
sources in our radio catalog, we expect to have a positional 
uncertainty of < 3 arcsec. 

The results of the Monte-Carlo simulations for the two 
source fitting packages are shown in Fig. 15. Plotted in 
linear-log scales is the fitted peak flux fractional difference 
dS/S = {Sinput — Sfitted)/S fitted, as a function of source 
strength. The two solid curves are quadratic least-square 
flts to the one-standard-deviation level of logj^Q(|dS'/5'j). The 
upper and lower bounds were fitted separately to better rep- 
resent the complex behavior of the source fitting algorithm 
in the presence of noise. The uncertainty has been set to a 
absolute minimum of 5% for strong sources. Overall, the es- 
timated uncertainty bounds > 90% of the data. The deficit 
of data points with low fitted flux densities {dS/S < 0) 
at faint source strengths is the well-known property of ra- 
dio source fitting in the presence of noise, as described in 
[Condon (1997). 

For the VSAD results in our present work, we have made 



^ The fitted source strength (us) is defined as the ratio of the 
fitted peak fiux density {Sp^ak) to the estimated local noise ((7„, 




0.6 



§2.2). 



Figure 14. The angular separation between the input and fitted 
source position as a function of fitted source strength, in linear-log 
scales. The solid curves indicates the fitted estimate of the posi- 
tional fitting uncertainty as a function of fitted source strength. 
The vertical dashed-lines denotes the various fitted signal-to-noise 
(c) level. 



the assumption that the fractional uncertainty in fitted peak 
flux density is representative of the fltted integrated flux 
density. 



APPENDIX II: ATCA AND NVSS FLUX SCALE 

Radio sources have an extremely wide range of luminosi- 
ties and usually requires a combination of surveys at var- 
ious sensitivity and aerial coverage to properly sample the 
radio population. In this section, we consider the possibil- 
ity of combining radio data from a survey with the ATCA 
and the NVSS from the VLA. We compared the source flux 
density measured in our survey and those in NVSS. With 
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Figure 15. The fitted peak flux deviation as a function of fit- 
ted source strength for the Monte-Carlo simulations, in linear- 
log scales. The solid curves indicates the fitted estimate of the 
fractional source fitting uncertainty as a function of fitted source 
strength. The vertical dashed-lines denotes the various fitted 
signal-to-noise (it) level. 



the on-line NVSS catalog browser^" we extracted the 157 
NVSS sources, brighter than 2.5 mjy, lying within our cur- 
rent ATCA-2dFGRS mosaic region. NVSS has a resolution 
of 45 arcscc and this is typically the minimum angular sep- 
aration between sources in the NVSS catalog. With higher 
angular resolution, some NVSS sources are resolved in our 
ATCA survey. To locate these 'decomposed' sources, we 
matched the 157 NVSS positions to our 5cr VSAD list using 
a search radius of 45 arcscc, centred at the NVSS source po- 
sitions. For matches with multiple components we represent 
the ATCA flux by simply adding the fitted integrated flux of 
the components. We used the NVSS flux uncertainty (pro- 
vided with the catalog) and our estimated flux uncertainty 
(derived with Montc-Carlo simulations, sec Appendix I) as a 
measure of the degree of deviation (the "sigma- level" ) of the 
data from the prefect one-to-one correlation. For instance, 
a source with NVSS flux Snvbs ± Anvss and correspond- 
ing ATCA flux 5atca ± Aatca, we defined the deviation 
"sigmarlevel" (Ao-), such that, for Snvss > "Satca, 

iSnvss — Act Anvss = Satca + Act Aatca, 

or in general, 

A _ Satca — Snvbs 
Anvss + Aatca 

We visually examined all NVSS positions with Act > 1.5 and 
those without an ATCA match. We found for a few extended 
sources, the Gaussian fitting method was inappropriate and 
we adopted the sum of the pixel values as the flux measure- 
ment. There are also a few cases where the source structure 
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Table 8. ATCA-NVSS cross identification results. 



Type 


Act: <1 


1-2 


2-3 


>3 


Total 




(1) 


72 


29 


7 





108 


(69%) 


(2) 


16 


3 








19 


(12%) 


(3) 


1 











1 


«1%) 


(4) 


7 


2 





1 


10 


(6%) 


(5) 


9 











9 


(6%) 


(6) 








6 


4 


10 


(6%) 


Total 


105 


34 


13 


5 


157 






(67%) 


(22%) 


(8%) 


(3%) 






(7) 





6 


1 


2 


9 





(1) Single ATCA component match. 

(2) Double ATCA components match. 

(3) Triple ATCA components match. 

(4) Extended structure requiring pixel summation. 

(5) Special cases. 

(6) Undetected in ATCA. 

(7) ATCA sources > 2.5 mJy without NVSS match (~ 6%). 




• Type(1) 

* Type (2) 
□ Type (3) 

1 : O Type (4) & (5) - 



1 10 100 

Figure 16. The NVSS and ATCA flux measurements of sources 

matched within our survey region. The difl'erent types of sym- 
bols denoted in the figure correspond to the matched-types listed 
in Table 8. The solid line is the least square fit to the data, 
the gashed-line is the one-to-one correlation between NVSS and 
ATCA fluxes. The dotted Unes notes the 2.5 mJy NVSS limit 
adopted. 



required some special care. The results of the cross identi- 
flcation process is listed in Table 8 and plotted in Fig. 16. 

From this small study, we did no observe any signs of signif- 
icance systematic differences between the two radio surveys, 
in particular, there do not seem to be any systematic lost of 
total radio flux density in the ATCA survey for intrinsically 
extended sources (i.e. sources resolved by ATCA). 



http: / /www.cv.nrao.edu/nvss/NVSSlist.shtml 



Table 3. Sample table of 2dFGRS-radio identifications and upper-limits. 



2dF position (J2000) 
oc 6 




2dF name 






^3 


Q 


class 


FM* 


VSAD fitted position (J2000) 


2dF-radio S^^^^ Sj^j^, 
Sep (") (mjy) (mjy) 


Flux u 
(%) 


ncertainty 
(mJy) 


(mJy) 




23^33" 


"'34.07" 


— 27° 44^41 


1" 


TGS192Z064 





0547 


17.27 


4 


a+e 


4 


- 


- 


- 


- 


0.8759 


- 


23"^ 33" 


136.31^ 


-28° ie/02 


2" 


TGS270Z051 





1042 


18.02 


4 


abs 


4 


- - 


- 


- 


- 


0.8026 


- 


23^33" 


137.29" 


-27°43'20 


1" 


TGS192Z058 





0546 


16.82 


5 


a+c? 


4 


- - 


- 


- 


- 


0.7719 


- 


23^33" 


'37.57= 


-28°00'42 


a" 


TGS270Z050 





1055 


17.20 


4 


a+e 


4 


- 


- 


- 


- 


0.7551 


- 


23^33" 


"'39.31" 


-28° 10'03 


3" 


TGS192Z060 





1073 


18.61 


4 


SF 


4 


- - 


- 


- 


- 


0.7094 


- 


23^33" 


"'39.96" 


-27°5l'06 


5" 


TGS192Z059 





0661 


18.94 


4 


SP? 


4 


- - 


- 


- 


- 


0.7028 


- 


23^33" 


142.84" 


-28° 16'04 


1" 


TGS192Z061 


n 


1042 


18.89 


4 


aba? 


4 


- - 


- 


- 


- 


0.6309 


- 


23^33" 


142.94" 


-28° 15'48 


2" 


TGS270Z048 





1042 


18.45 


4 


abs? 


4 


- 


- 


- 


- 


0.6260 


- 


23^33" 


'43.63" 


-28°2l'08 


s" 


TGS192Z062 





1054 


18.73 


4 


a+e 


4 


- - 


- 


- 


- 


0.6702 


- 


23*^33" 


"'48.94" 


-27° 25^25 


1" 


TGS192Z055 





1504 


19.31 


4 


em 


4 


- 


- 


- 


- 


0.5655 


- 


23^33" 


"'51.17" 


-27°53'20 


8" 


TGS192Z056 





2159 


19.38 


1 


??? 


4 


- - 


- 


- 


- 


0.4698 


- 


23^33" 


152. 60" 


-28°25'05 


2" 


TGS270Z042 





1064 


18.82 


4 


a+c 


4 


- - 


- 


- 


- 


0.5768 


- 


23^33" 


152.97" 


-27°39'41 


q" 


TGS192Z054 





1049 


17.46 


5 


abs 


4 


- - 


- 


- 


- 


0.4482 


- 


23^33" 


154.13" 


-28°20'52 


q" 


TGS192Z057 





2162 


18.69 


1 


??? 


4 


- - 


- 


- 


- 


0.4767 


- 


23*^33™58.07^ 


— 27° 42'07 


0" 


TGS192Z052 





2053 


19.35 


1 


??? 


4 


- - 


- 


- 


- 


0.3877 


- 


23 33'' 




-27° 23^23 


q" 


TGS192Z051 





1012 


18.14 


4 


abs 


4 


- 


- 


- 


- 


0.4536 


- 


23^ 34" 


"'01 04" 


-27°42'59 


5" 


TGS192Z053 





2157 


18.24 


1 


??? 


5 


- - 


9.500 54.9 


8 


4.640 


0.3581 


27 


23^ 34P 


"■Ol 28" 


-28°25'53 


g" 


TGS270Z037 





1537 


18.94 


4 


abs 


4 


_ _ 




- 


- 


0.4709 


- 


23^ 34" 


''01 40" 


-27° 17^47 


3" 


TGS192Z049 





1057 


18.85 


4 


em 


4 


_ _ 




- 


- 


0.5139 


- 


23^ 34P 


''01 43" 


-28°33'08 


1" 


TGS270Z039 





1049 


18.48 


5 


abs 


4 


- 


- 


- 


- 


0.7917 


- 




''01 90" 


-28°27'39 


9" 


TGS270Z038 





1045 


19.05 


5 


SF 


4 


- - 


- 


- 


- 


0.5096 


- 


23'>34= 


102 61" 


-27°36'31 


1" 


TGS192Z050 





1522 


19.10 


4 


abs 


4 


- - 


- 


- 


- 


0.3436 


- 


23'>34= 


'03 35" 


-28° ir/46 


3" 


TGS270Z036 





2764 


19.38 


4 


abs? 


3 


- - 


0.459 


157 


0.721 


0.3447 


1 


2 h 3^11 


'04 46" 


-28°23'l9 


5" 


TGS270Z032 





1500 


18.95 


4 


??? 


4 


- 


- 


- 


- 


0.3941 


- 


23l>34= 


104 77" 


-28°24'55 


8" 


TGS270Z033 





1222 


18,78 


5 


abs 


4 


- - 


- 


- 


- 


0.4142 


- 


23'>34= 


'07 88" 


-28°35'41 


6" 


TGS271Z031 





1539 


18.72 


4 


em? 


4 


- - 


- 


- 


- 


0.8620 


- 


23'>34= 


'08 24" 


-27°22'26 


0" 


TGS192Z047 





1002 


19.36 


4 


SF? 


4 


- 


- 


- 


- 


0.3667 


- 


2 h 3^11 


'09 80" 


-27° 18^32 


q" 


TGS192Z045 





1006 


19.17 


4 


SF 


4 


- - 


- 


- 


- 


0.4066 


- 


23'>34= 


'09 89" 


-28°30'52 


3" 


TGS270Z031 





0863 


18.74 


5 


abs 


3 


- 


0.565 


166 


0.935 


0.5142 


1 


23^34" 


1 12.46" 


-27° lO'lS 


8" 


TGS192Z185 





1203 


19.08 


4 


a+c? 


4 


- - 


- 


- 


- 


0.6626 


- 


2311341] 


1 15.03" 


-28°26'34 


n" 


TGS270Z030 





1219 


18.88 


5 


SF? 


4 


- - 


- 


- 


- 


0.3558 


- 






-27°28'37 


1" 


TGS192Z046 





1013 


17.92 


4 


a+e 


4 


- - 


- 


- 


- 


0.2770 


- 




'15 74" 


-27°27'l2 


7" 


TGS192Z042 





1947 


19.28 


1 


??? 


4 


- - 


- 


- 


- 


0.2810 


- 


23113411 


'17 09" 


~27°53'32 


q" 


TGS192Z044 





1598 


19.00 


4 


??? 


4 


- - 


- 


- 


- 


0.2465 


- 


23^34" 


'20.88" 


-28°32'31 


2" 


TGS270Z029 





1051 


18.77 


4 


SF 


4 


- 


- 


- 


- 


0.4563 


- 


23^34" 


'20.91" 


-28°09'43 


\" 


TGS270Z028 





2159 


19.03 


1 


??? 


4 


- - 


- 


- 


- 


0.2301 


- 


23^34" 


'23.46" 


-27° 15'57 


o" 


TGS192Z040 





0882 


18.06 


4 


abs 


4 


_ _ 




- 


- 


0.3524 


- 




^23 60^ 


-27°38'06 


6" 


TGS192Z041 





1053 


19.27 


4 


abs? 


4 


_ _ 




- 


- 


0.2183 


- 


23^34" 


^23.75^ 


-27°58'l0 


2" 


TGS270Z026 





2120 


19.02 


4 


a+c? 


4 


- 


- 


- 


- 


0.2149 


- 


23^34" 


"^26.96^ 


-28°05'25 


7" 


TGS270Z027 





2159 


17.84 


1 


? ? ? 


4 


- - 


- 


- 


- 


0.2030 


- 


23^3411 


^32.20^ 


-27°55'l0 


0" 


TGS192Z037 




- 


18.64 


1 


??? 


4 


- - 


- 


- 


- 


0.1836 


- 


23^34" 


■^32.30^ 


-27°3l'l3 


1" 


TGS192Z036 





0860 


19.00 


4 


abs? 


4 


- - 


- 


- 


- 


0.1957 


- 


23^34" 


■^33.59^ 


-28°40'l9 


7" 


TGS270Z025 





1041 


17.82 


5 


abs 


4 


- 


- 


- 


- 


0.8460 


- 


23^34" 


^35.18^ 


-27°59'43 


0" 


TGS192Z038 





0866 


17.25 


4 


abs 


4 


- - 


- 


- 


- 


0.1743 


- 


23I13411 


"^38.77^ 


-28° 17'l7 


1" 


TGS270Z022 





1836 


19.20 


5 


abs 


4 


- - 


- 


- 


- 


0.1797 


- 


23^3411 


^39.41^ 


-28°0l'37 


3" 


TGS192Z034 





2121 


19.34 


4 


SF? 


4 


- - 


- 


- 


- 


0.1619 


- 


23^34" 


"40.05^ 


-27° 49^47 


9" 


TGS192Z033 





0278 


16.87 


4 


SF? 


4 


- 


- 


- 


- 


0.1600 


- 


23^^34" 


■^41.24^ 


-28°3l'37 


1" 


TGS270Z023 





0370 


19.07 


4 


??? 


4 


- 


- 


- 


- 


0.3015 


- 


23^34" 


^45.31^ 


-28°34'52 


7" 


TGS271Z028 





1054 


18.12 


4 


SF? 


3 


- - 


0.446 


161 


0.719 


0.3641 


1 


23^3411 


^45.98-^ 


-27°51'34 


0" 


TGS192Z032 





2161 


19.31 


1 


??? 


3 


- - 


0.166 


164 


0.272 


0.1444 


1 


23^34" 


^49.29^ 


— 28° 13 26 


9" 


TGS192Z031 





2159 


18.21 


1 


??? 


4 


- - 


- 


- 


- 


0.1470 


- 


23^34" 


■^50. 81^ 


-28°38'32 


4" 


TGS270Z019 





1040 


18.93 


3 


em? 


2 


- - 


2.532 


52 


1.328 


0.4878 


5 


23^34" 


■^51.62^ 


-28°4l'20 


6" 


TGS270Z020 





1503 


17.94 


4 


abs 


4 


- - 


- 


- 


- 


0.7119 


- 


23^3411 


^55.19^ 


-27°0l'47 




TGS192Z166 










??? 












8262 




23^34" 




-28°38'l8 


7" 


TGS270Z018 





1497 


18.65 


4 


abs 


4 


- - 


- 


- 


- 


0.4269 


- 


23''35'' 


^01.28^ 


-28°43'l3 


h" 


TGS271Z026 





1494 


19.36 


4 


abs 


4 










0.8413 




23'*35" 


■^01.81^ 


— 28°43'02 


3" 


TGS271Z025 





1500 


18.44 


5 


abs 


4 










0.8085 




23''35'' 


■^03.43^ 


-27° 18'50 


2" 


TGS193Z018 





1260 


19.34 


4 


SF 


3 




0.199 


164 


0.327 


0.1730 


1 


23^35" 


^04.05^ 


-27° 17'l6 


a" 


TGS192Z029 





1273 


19.24 


4 


abs 


4 










0.1850 




23I135" 


"^04.65'' 


-27°25'l5 




TGS192Z030 





1246 


17.76 


1 


??? 


4 










0.1350 




23''35"' 


^06.53^ 


-27° 14'47 


q" 


TGS192Z028 





0510 


17.76 


5 


SF 


4 










0.2055 




* The 


fitting n 


ncthod used tc 


dcri 


/c the radio flux 




;hc upp 


3r-limit, 


(1) 


VSAD 4(7 n 


latch, see 


§3.2; (2) pixel sum, see §3.3; (3) 


non-identification upper-limit, se 


e §3.4; 









(4) IMFU upper-limit, see §3.5; (5) local noise upper-limit, see §3.5 
The estimated VSAD fitted positional uncertainty in arcsce, see Appendix I. 
The estimated local noise, see §3.1.2. 
The signal-to-noise ratio (Spg^j^ /cts ). 



